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Abstract 

We supplement Maximum Likelihood methods with a Monte-Carlo sim- 
ulation to re-investigate the SN1987A neutrino burst detection by the 1MB 
and Kamiokande experiments. The detector simulations include background 
in the the latter and "dead-time" in the former. We consider simple neu- 
trinosphere cooling models, explored previously in the literature, to explore 
the case for or against neutrino vacuum mixing and massive neutrinos. In 
the former case, involving kinematically irrelevant masses, we find that the 
full range of vacuum mixing angles, < sin^ 29v < 1, is permitted, and the 
Maximum Likelihood mixing angle is sin^ 29v = .45. In the latter case we 
find that the inclusion of "dead-time" reduces previous m^^ upper bounds by 
10%, and supplementing the Maximum Likelihood analysis with a Monte- 
Carlo goodness-of-fit test results in a further 15% reduction in the m^^ upper 
limit. Our 95% C.L. upper limit for m^^ is 19.6eV, while the best fit value is 
~ OeV. 
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Introduction 

Galactic neutrino astronomy began in 1987 with the observation of 20 
neutrinos from the supernovae burst SN1987A in the Large Magellanic Cloud 
(LMC). Two terrestial detectors, 1MB [0 and Kamiokande 0, found un- 
equivocal evidence for supernovae neutrino events with the former collab- 
oration claiming detection of 8 SN1987A events, and the latter 12. This 
momentous observation generated enormous excitement in the scientific com- 
munity, and of course a plethora of papers soon followed. We have returned 
to this subject for three reasons. First, apparently only the 1MB collabora- 
tion took into account the "dead-time" in the 1MB detector when comparing 
these observations with theory. Also we believe that a more or less model- 
independent approach to the question of vacuum mixing in SN1987A should 
be done (and with more rigor than in [^). Finally, we have in hand a detailed 
Monte Carlo code for generating the predicted signal in light water detectors 
which was created for the purpose of exploring the nature of a galactic neu- 
trino signal 0] , but which can also be used to accurately model the signal for 
SN1987A. Using this code we felt that it might be possible to improve upon 
the neutrino mass limits derived previously for the SN1987A neutrino burst. 

Before describing our analysis, it is worthwhile briefly reviewing the most 
recent results on neutrino mass and mixing constraints from SN1987A in or- 
der to comment on the improvements incorporated in our present analysis. 
To date, the most comprehensive statistical analysis of the signal was per- 
formed by Lamb and Loredo [^j , who used a Maximum Likelihood technique, 
and were the first to incorporate the background event rate in a likelihood 
function in order to account for at least one of the Kamiokande events which 
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was probably a background event (which several previous analyses had sim- 
ply discarded). To facilitate comparison with their results (in the case of a 
significant neutrino mass) we exploited essentially the same formalism, but 
made several minor additions. Primarily, we included the fact that the 1MB 
detector had a significant amount of "dead-time" (~ 13% over theburst |jl|) 
following interactions in the detector by cosmic ray muons or SN1987A neu- 
trinos. The reason for the large dead-time, 35ms/interaction, is due to the 
data-acquisition software This 35ms dead-time would of course not be 
problematic for the original purpose of the 1MB detector, measuring proton 
decay. 

A few investigators p, 0] allowed independent offset times (the time-delay 
between the first SN1987A neutrino arriving at the earth and the first detec- 
tion in either Kamiokande or 1MB) for each detector. We noticed that in the 
case of a significant neutrino mass, when the offset times, toff, are impor- 
tant, the Maximum Likelihood value t^j™ (Kamiokande offset time) at the 
m,^^ upper hmit of Lamb and Loredo seemed unacceptably long. This sug- 
gested extending the Maximum Likelihood approach. The Maximum Like- 
lihood method does not test a model, but rather tests the allowed range of 
parameters given a specific model. Thus should the Maximum Likelihood 
method locate a parameter value which may seem otherwise unlikely one 
must use other methods, such as a Monte-Carlo, to test the model. Here 
we were able to exploit the power of the comprehensive light water neutrino 
detector Monte Carlo code previously written to examine various features of 
possible future galactic neutrino burst signals 0]. This code incorporates all 
aspects of the detector in order to generate a realistic signal, given a specific 
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supernova burst model. Using this Monte Carlo, we can show, as we describe 
in more detail below, that the Lamb and Loredo Maximum Likelihood offset 
time for Kamiokande of 3.9 sec would be expected to occur in less than 1 
out of 400 cases given the other Maximum Likelihood SN1987A parameters, 
such as binding energy, emission timescale, etc. 

We next turn to the issue of neutrino mixing constraints. There has 
recently been a model-dependent derivation of bounds on vacuum oscillations 
from the SN1987A data [Q, with 5 specific explosion models considered. All 
models fit the no-mixing scenario and, perhaps not surprisingly, when the 
amount of neutrino mixing was increased a Kolmogorov-Smirnov test led to 
a lower bound of sin^26'y < .7 — .9. The upper bound of .7 would exclude 



the "just-so" solution to the solar neutrino problem [jTO| and much of 
the large angle region of the MSW solution. The severity of this bound is 
surprising, given the sparsity of the observed signal, so we decided to examine 
this issue in some more detail. The authors of ^ recognized the fact that 
their result was model dependent, but just how model dependent was not 
clear. The ability to explore the neutrino signal with our Monte Carlo makes 
it very easy to sample supernova model space. We will show below that with 
a minimally model-dependent approach, maximal vacuum mixing actually 
fits the data better (greater likelihood) than no mixing. 

Finally, we note that the neutrino mass limit we derive here has already 



been superceded by direct laboratory probes |]rT|. Nevertheless, the utility 
of exploiting a galactic supernova burst to constrain neutrino masses and 
mixings remains of great interest, and the techniques we examine here thus 
remain important to explore. Namely, SN1987A remains, even 7 years later. 



4 



an important, and unique test case if we are to attempt to fully exploit the 
information which may be available in the next observed supernova neutrino 
burst. 

The Minimal Model 

We exploit here two "minimal" models. The version we use for bounding 
the neutrino mass and the one used for bounding the vacuum mixing angle 
differ in that the latter has an extra time constant in the neutrino spectra, 
which we will discuss later. In the former case we follow Lamb and Loredo 
and assume a simple exponential cooling model. In this case the supernovae 
is characterized by a Fermi-Dirac neutrino spectrum and 3 parameters; a 
maximum initial temperature T^^, a cooling time-scale Tc, and a, related to 
the size of the neutrinosphere. 



With it!io the radius of the neutrinosphere in units of 10 kilometers and Dr,o 
the distance to the LMC in units of 50 kiloparsecs. Alternatively one can 
view a as a relation for the supernovae binding energy with the additional 
assumption that there is an equipartion of the binding energy carried away 
among the 3 flavor states times 2 spin states of the emitted neutrinos. 



with £^53 the neutron star binding energy in units of lO^^ergs and Tp^ given 
in the simple exponential cooling model by 
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exp(-t/4rc). 



(3) 



The model is also characterized by the neutrino mass, rriy^, and an ad- 
ditional parameter for each detector, toff., the offset time. The offset time 
is particularly important for the massive neutrino case where the neutrino 
mass causes a delay-time in the arrival of the neutrinos: 



with niy^ in units of eV and EueVi the incident neutrino energy, in MeV. 
Independent offset times are needed for each detector 0] because of a problem 
with the Kamiokande clock during the time that the SN1987A neutrino burst 
passed the earth. The offset times play a major role in constraining non- 
zero neutrino masses, due to the difficulty, when maximizing the likelihood 
function for the Kamiokande detector, of reconciling a few early low energy 
neutrino events, which would then imply a large offset time for non-zero 
mass, with following high energy events, which would favor a small offset 



The second version of our minimal model, in the case of vacuum mixing 
and nearly massless neutrinos, introduces an additional parameter for the 
anti- neutrino temperature, Tp^. As far as kinematics are concerned, in this 
part of the analysis we assume effectively massless neutrinos (the time delays 
introduced by the very small masses of interest in this case are irrelevant), and 
we assume two state mixing, v^, and v^. Vacuum mixing implies the neutrino 
spectrum at the earth is a mixture of the 2 original spectra according to 




(4) 



time p|. 
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(fN'^^ydEdt = (1 - .5 sin^ 2ev)d^N''^/dEdt + .5 sin^ 2evd^N^^ /dEdt (5) 



e 



Qualitative arguments suggest that T^^ ~ 2T^^ due to the fact that the 
have additional interactions from neutrons and charged current interactions, 
while z/^ has only neutral-current interactions in the supernovae environment. 
Thus the z/^ are emitted from deeper in the star and hence their spectrum is 
characterized by a hotter temperature. The temperature evolution of the Ue 
is potentially slightly more complicated: initially before neutronization of the 
star Tp^ ~ T^^ is expected. However as the star neutronizes and the reaction 
i'e + P ^ n + becomes rare then Tp^ is expected Our model 



parameterizes this phenomenon by introducing an additional time-constant 
which smoothly and symmetrically takes Tp^ from T^^ to Ty^ = 2Ty^ in time 
2r2. 



Tp^ = {1. + .5tanh [{t - t^)^ / r^)]} T.^t) (6) 

This function is constructed so that when t advances to T2 we have Tp^ = 
.5(Tj^^ + Ty^). Introducing the extra parameter, r2, allows us to take a con- 
servative approach to the consideration of vacuum oscillations. Also we take 
Tj^g and Ty^ to have the form of eq.^ (such that if the Maximum Likelihood 
T~2 ^ Tc this will reduce Tp^ to the form used for massive neutrinos). 

Another point worth mentioning here is that we continue to partition the 
binding energy equally among all species here so that the factor a, which 
sets the scale of the fluxes. 



d^N/dEdt oc 
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is different for tlie 3 species, and i/^. Tliese are related by, 

a.^^oc.^iijTtdt/ jrtM) (7) 

Witll Ui = Uejl^f,. 

As we liave indicated, an improvement in our Maximum Likelihood method 
compared to the Lamb and Loredo work is the inclusion of dead-time for the 
1MB detector. The reason one expects this may have an effect is due to 
the fact that the Kamiokande data favors a cooler, less energetic supernovae 
than does the 1MB data. Our Monte-Carlo work indicates that if one uses 
the 1MB data to locate a set of SN1987A parameters, these same parame- 
ters typically predict many more events in Kamiokande ~ 20 , with a much 
higher average energy ~ 22MeV, than were seen ( 12 events and 15 MeV, 
respectively) . Thus one of the reasons why the Maximum Likelihood analysis 
can provide a reasonably localized parameter space for the combination of 
the 1MB and Kamiokande data is the tension between the fits for the two 
separate data sets. Since including the dead-time in 1MB will favor an even 
more energetic supernovae it should exacerbate the existing tension leading 
to stronger constraints on parameters derived from the Maximum Likelihood 
method. One of the purposes of this letter is to show how much the 13% 
dead-time changes our conclusions from those of Lamb and Loredo . 

The dead-time, ta — 35ms, was handled in different ways in the Monte- 
Carlo and Maximum Likelihood methods. In the prior case it is straightfor- 
ward. We use a Poisson distributed random number generator to simulate 
the known 2.7Hz muon event rate (this requires a 3Hz incident rate since 
dead-time affects this measurement as well) starting at t — Isec. With the 
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Monte-Carlo neutrino events and muon events in temporal order we then re- 
move any neutrino events occurring within 35ms of a previously "detected" 
muon or neutrino. 

For the Maximum Likelihood method, we modify the spectral rate ac- 
cording to, 

d^N/dEdt ^ (1 - Pd{t - td, t))d^N^/dEdt (8) 

where Pd{t — t^i, t) is the probability that either a muon event or a neutrino 
event ocurred in the interval from t—td to t and d'^N^/ dEdt is the spectral rate 
without deadtime. The probability, P^, of an interaction which causes dead- 
time is decomposed as follows Pd = Pd^Pd„, using the poisson probability 
that there were events from t — t, assuming a rate F. 

Pd = l- exp i-TSt) (9) 

For the muons F^ = 3Hz, 6t = 35ms. We approximate the neutrino induced 
dead time probability by evaluating the zeroth order contribution of the 
neutrinos to Pd at t — td/2, approximating d'^N^ /dEdt as constant during 
that short interval. Since td <^ Tc this is a good approximation. Thus, 
Fjy oc dN^^ and 6t = min{t, td). 

The data in our Maximum Likelihood code not specifically mentioned 
above, such as fiducial detector volumes, the parameters of the detector 
resolution functions, the energies and times of the background events in 
Kamiokande, follow the treatment in Lamb and Loredo , and are not re- 
peated here. We use the standard likelihood function, see for example 0, |]. 
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The Monte-Carlo program is a modified version of tlie one described in [^]. 
Tlie parameters wliicli describe tfie detector efficiencies, resolution functions, 
the form of the neutrino temperature etc, have been set to be identical to the 
ones for the Maximum Likelihood analysis. This code (originally designed 
for 0(1000) events) is fairly sophisticated and includes the interactions of 
neutrinos other than i/g, and neutrino scattering from oxygen nuclei in the 
detector. The additional types of interactions in this code, and the more 
careful treatment of the dominant reaction v^P ne~^, which includes nu- 
cleon recoil effects, results in a small increase, ~ 3%, in < > compared to 
the Maximum Likelihood code estimate. The difference is insignificant for 
the SN1987A events however, as will become obvious. 
Analysis of Results 
(a) Massive Neutrinos 

We first consider the limits on massive neutrinos. The initial step is 
to find the best fit Maximum Likelihood parameters for a, T^^, Tc, toj™ and 
as a function of m,^^. The effect of this procedure is to project the log 
likelihood onto the m,^^ axis. From this projection, shown in Figure 1, we 
can find the 95% confidence limit from, 

ln£^--.5xL/.05), (10) 

where we have 6 degrees of freedom [dofj for the chi-squared distribution 
in the present instance. We will denote the Maximum Likelihood value of 
a parameter by adding the subscript L+, and the value of a parameter at 
the 95% confidence boundary with some or all (which will be clear from the 
context) of the other parameters at their Maximum Likelihood values by 



10 



adding the subscript L . In this notation, from Figure 1, 



{m,jL- = 23eV. 



(11) 
(12) 



The hkehhood function is extremely flat below m,^^ = 2eV so this result 
does not strongly favor an identically zero neutrino mass. Our value for 

{my^)i- is 8% lower than the Lamb and Loredo result {m^Ji^- = 25eV, the 
entire difference being due to the dead-time correction in our 1MB likelihood 
function. 

There is an additional constraint we may use in the analysis. Con- 
sider Figure 2 wherein (t^ff) l+ and (fjff) l+ are shown as a function of m,^^ 
(with (tc)l+ included for comparison). Note that {tgfJ')L+ reaches 4.2s at 
rrii,^ = 23ev (where (^o™/) l+ = ^s) . The Maximum Likelihood offset time for 
Kamiokande seems extraordinarily long, especially in light of the fact that 
{'toff)L+ exceeds (tc)l+ for m,,^ > 21.7 eV. 

To test our intuition in this regard, and to discover the acceptable range 
of f^ff one can use our Monte-Carlo code to find, given {0)1+, {T^Jl+ and 
{'Tc)l+, for a particular m„^, the probability, P{toff < t\N^^,m,y^), that the 
offset time in a particular detector will not exceed a certain value. 

Note that we are interested in values of the mass parameter in the range 
{^ue)L+{OGV) < < i^ue) L- {23eV) . Wc thus first determine, using the 
Maximum Likehhood method, {t'^'^J') l- ■, the minimum acceptable t'^jj', sub- 
ject to the constraint with all the other parameters free. This is dis- 
played in Figure 3. Then we use our Monte-Carlo to construct P{tgfY < 
t\N'"'"',m^J. If we find that {t';'}f )L- is ruled out at the 95% CL by this 
probability distribution, this then implies, since the likelihood function rules 
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out any smaller offset time, that the m^^ corresponding to this value is at 
least as unacceptable at this level. 

We finally turn to the constuction of P{t'^'}f < t|A^*^""*, m^J. Recall that 
in our detector simulation we include dead-time for 1MB and background 
for Kamiokande. Also note that the time of the first event depends (more 
strongly as the number of events is decreased) on the number of events de- 
tected, and that this number is not fixed by our Monte-Carlo code, which 
temporally simulates the neutrino burst and detection. Therefore we re- 
quire Monte-Carlo runs which result in the desired number of events for 
each detector. Then we rank the times of the first event of each such run 
and generate a cumulative probability distribution for the time of the first 
event. To improve the statistics, while conserving computer time, we choose 
a range, N = N'^'^^ ± 1, about the desired number of events. (For this pur- 
pose N^"-"^ — 16, including background, and N^'^^ — 8.) ( Both detectors, 
for the parameters of interest here, have a fiat distribution for the expected 
time of the first event in the neighborhood of the number of events actually 
observed. ) We use 1000 Monte-Carlo runs of each detector to acquire the 
data for the construction of Pit^fj < t\N'^^\ m^J. TypicaUy about 20 - 30% 
of the Monte-Carlo runs fall in the accepted range of N'^^. 

In Figure 4a we plot P(iJ*^ < t\N 2leV) . Also shown are {tlff)L+, 
and (47/)l+> for m^^ = 21eV . Note that while (47/)l+ = -Qs is located near 
the mean of the distribution, {t^'ff)L+ — 3.8s is in the tail. In Figure 4b 
P(tf//" < t|A^^«"*, 19.6ey) is plotted with {t^oW)L+ and {t^o7f)L- indicated. 
Since {t';Jf)L- > t*, where t* is defined by P{t';'}f < t*|iV^'"™, 19.6eV) = .95, 
^i/, > 19.6eV is excluded at the 95% confidence level. 



12 



We thus find an additional decrease of 15% in the m,y^ upper hmit deriv- 
able using the Monte-Carlo generated probabihty P{t^off < t\N'^°''^,m^J in 
addition to the Maximum Likelihood procedure. This is a significant factor, 
and further underscores the utility of Monte Carlo simulation of the data. 
When combined with the effect of incorporating deadtime in the 1MB de- 
tector, which fortuitously plays a significant role because of the paucity of 
observed events in Kamiokande, we have been able to reduce the upper limit 
on the m,^^ mass by over 25% compared to previous analyses. 

(b) Vacuum Mixing and Nearly Massless Neutrinos 

Next we turn to our results for vacuum mixing. In this case the Maximum 
Likelihood values for the offset times are for all values of sin^ 2^y, therefore 
we do not have to Monte-Carlo the neutrino burst. In Figure 5 ln£(sin^ 2^y) 
is displayed for the range < sin^ 29v < 1. (We ran the Maximum Likelihood 
code with and without including the offset times. The difference in In C never 
exceeded .02%. Thus the offset times are irrelevant parameters, and we 
consider only 5 dof.) In Figure 5 the likelihood function peaks at sin^ 29y = 
.45. However, the likelihood function is relatively fiat over the entire range 
so non-zero mixing is only marginally preferred. The likelihood ratio for 
sin^ 29v = 0.45 compared to sin^ 29v = is 5.5. 

Also of interest, perhaps to supernovae model builders, is our Maximum 
Likelihood extraction of neutrinosphere temperatures. In Figure 6a-c we 
display Tu^,Tp^ and T^^ for sin'^ 2^^^ = 0, .45 and 1 respectively. The main 
feature is that in all cases Tp^ T^^ gradually, with T2 = 9.13, 8.72 and 8.35 
respectively. This long timescale is something of a surprise. 

A final question is whether the admission of constraints on SN1987A pa- 
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rameters, other than those purely obtainable from the neutrino data alone, 
would allow one to further limit sin^ 26'y. In Figure 7a-c we show: the Max- 
imum Likelihood neutron star binding energy, E^^, in units of 10^^ ergs; the 
intitial electron neutrino temperature, T°^, in MeV; and the cooling time- 
scale, Tc, in seconds. The entire range of the latter seems acceptable based 
on estimates from supernova models. One may ask whether constraints such 
as E^^ < 4.5 or T^^ > 3MeV, would limit sin^26'y. From Figure 7 it appears 
that this could be the case. 

To address this question we find the 95% confidence level regions in the 
E^^ , plane for several values of sin^26'y. The remaining parameters are 
all permitted to find their Maximum Likelihood value^. In Figure 8 we 
display the results for 6 values of sin^ 29v from .1 to 1. The outer contour in 
each box is the 95% confidence limit (5 degrees of freedom), while the inner 
contours; 50%, 25%, and 10% CL. , are shown to allow the reader to assess 
the character of the surface. It is apparent from Figure 8 that no reasonable 
E^2, provides a further solid constraint on sin^2^v'- 

We have constructed Table 1 to present the maximum allowed sin^ 29v 
at 95% confidence level in terms of a given minimum permissible T^j^, which 
we designate Tf; this is understood to refer to a limitation on the electron 
neutrino temperature provided independently of the neutrino data, such as 
may arise from supernovae modeling. If one cannnot bound from below 
then the parameter sin^ 29v cannot be constrained. In order to rule out large 
vacuum angle solutions to the solar neutrino problem (.7 < sin^2^y < .9 
[|]), using the SN1987A neutrino data, a rigorous argument that supernovae 

^In our formalism E^^ is determined by the combination of a, and Tc, thus in the 
present context we obtain a from the fixed values of E^^ and and the free value for Tc- 
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dynamics require > 4 MeV appears to be needed. The parameters 
in Table 1 are well within the typical range of 3-5 MeV in the supernovae 



model literature |J3|, |14 
Conclusions 

Our results demonstrate several important lessons for statistical analyses 
of constraints on neutrino properties from a nearby supernova neutrino burst, 
as well as refining these constraints for the observed burst from SN1987A. 
In the first place, while a Maximum Likelihood procedure can provide very 
powerful constraints on model parameters, it alone cannot address the ques- 
tion of whether these model parameters can be realistically achieved. When 
these parameters have to do with features of the observed burst, and not 
internal features of an underlying supernova model, then a Monte-Carlo pro- 
cedure such as we have devised 0] can prove to be very useful in further 
strengthening constraints on neutrino properties. 

Next, we have seen that the ability of the SN1987A burst signals in 
Kamiokande and 1MB to constrain a non-zero electron neutrino mass is in 
some sense fortuitous, due to the "tension" of the Kamiokande and 1MB 
data — in particular the apparent paucity of events in Kamiokande relative to 
1MB. For this reason, when we included deadtime in 1MB we were able to fur- 
ther extend the lever-arm in constraining m,y^. Our final result, m,y^ < 19.6eV 
is approximately 25% stronger than the previous best limit. 

Finally, we find that the ability of the combined SN1987A neutrino bursts 
to constrain neutrino masses does not at present extend to an ability to con- 
strain neutrino mixing angles in any model independent way. In particular, 
because of the uncertainty in the timescale for neutronization, without intro- 
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ducing strong model dependence — in particular contraints on E^^ and Tp^ — 
one cannot limit sin^ 29v from the SN1987A neutrino events. This argues 
against the claim made in [Q. It will be interesting to determine just how 
strong the constraints might become for a galactic supernova burst, and this 
issue is currently under investigation. 

We thank Steve Dye, Robert Svoboda and Martin White for useful con- 
versations. We also thank the 1MB collaboration for providing unpublished 
data. 



16 



References 

[1] C.B.Bratton et al. , Phys. Rev. D 37, 3361 (1988) 
R.M.Bionta et al. , Phys. Rev. Lett. 58, 1494 (1987) 

[2] K.S.Hirata et al. , Phys. Rev. D 38, 448 (1988) 
K.S.Hirata et al. , Phys. Rev. Lett. 58, 1490 (1987) 

[3] L.M.Krauss, Nature 329, 689 (1987) 

[4] L.M.Krauss, P.RomaneUi, D.Schramm, and R.Lehrer, NucL Phys. B. 
380, 507 (1992) 

[5] T.J.Loredo, and D.Q.Lamb, Ann. N.Y. Acad. Sci. 571, 601 (1989) 
[6] Steve Dye, private communication 

[7] L.F.Abbot, A.De Rujula, and T.P.Walker, Nucl. Phys. B. 299, 734 

(1988) 

[8] A.Y.Smirnov, D.N.Spergel, and J.N.Bahcall, preprint |hep-ph 9305204 
(1993) 

[9] P.LKrastev and S.T.Petcov, Phys. Rev. Lett. 72, 1960 1994 
P.J.Kernan, Ph.D. Thesis Ohio State University (1993) 
S.L.Glashow and L.M.Krauss, Phys.Lett. B190 199 (1987) 
V.Barger, K.Whisnant and R.J.N.Philhps, Phys. Rev. D 24 538 (1981) 
J.N.Bahcall and S.C.Frautschi, Phys. Lett. B29, 263 (1969) 
B.Pontecorvo, Sov.JETP, 26 984 (1968) 



17 



[10] An excellent review is found in J.N.Bahcall, Neutrino Astrophysics 
(Cambridge University Press, Cambridge, 1989). 

[11] Particle Data Group, Phys.Rev.BbO 1173 (1994) 

[12] E.W.Kolb, A.J.Stebbins, and M.S.Turner, Phys. Rev. D. 35, 3598 (1987) 
[13] R.Mayle, J.R.Wilson, and D.N.Schramm, Ap.J. 318, 288 (1987) 
[14] A.Burrows, Ann.Rev.Nucl.Sci. 40, 181 (1990) 



18 



Figure Captions 



Figure 1: The projection of the log Maximum Likehhood onto the tjIi,^ axis. 

Figure 2: The Maximum Likehhood offset times for the 1MB and Kamiokande 

detectors as mj,^ is varied. For comparison the Maximum Likehhood super- 
novae coohng timescale is also shown. 

Figure 3: The log likelihood as a function of the offset time in the Kamiokande 
detector for several values of mi^^. The horizontal line indicates the 95% C.L. 
boundary. 

Figure 4: Shown are Monte Carlo generated cumulative probability distri- 
butions for the time of the first event in a detector given the supernovae 
model parameters. In (a) the neutrino mass is 21 eV, and the distributions 
for 1MB and Kamiokande are shown. The Maximum Likelihood offset times 
are also indicated. In (b) the neutrino mass is 19.6 eV. The distribution for 
Kamiokande is shown, as are the Maximum Likelihood and 95% C.L. offset 
times. The short horizontal line is at P=95%. 

Figure 5: The projection of the log Maximum Likelihood onto the sin^ 29v 
axis. 

Figure 6: Temporal profiles of the Maximum Likelihood electron, anti-electron, 
and muon neutrinosphere temperatures for mixing angles of sin^ 29v = (a), 
sm'^29v = .45 (b) and sin^ 2^y = 1 (c). 
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Figure 7: As a function of neutrino mixing angle, the Maximum Likelihood 
neutron star binding energy in units of 10^^ ergs (a), initial electron neutri- 
nosphere temperature in MeV (b), and supernovae neutrinosphere cooling 
timescale in seconds (c). 

Figure 8: Maximum Likelihood Projections into the plane of the neutron star 
binding energy and initial electron neutrinosphere temperature for several 
choices of the neutrino mixing angle. The binding energy is in units of 10^^ 

ergs and the temperature in MeV. The contours are displayed at the 95%, 
50%, 25% and 10% Confidence Levels. For sin^ 2^y = 1 the 10% C.L. contour 
does not exist. For sin^ 29v — .1 the 25% and 10% C.L. contours do not exist. 
(See Figure 5). 

Table 1 : Maximum on the 95% C.L. boundaries of Figure 8 
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